Semiclassical Study on Tunneling Processes via Complex-Domain Chaos 
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We investigate the semiclassical mechanism of tunneling process in non-integrable systems. The 
significant role of complex-phase-space chaos in the description of the tunneling process is elucidated 
by studying a simple scattering map model. Behaviors of tunneling orbits are encoded into symbolic 
sequences based on the structure of complex homoclinic tanglement. By means of the symbolic 
coding, the phase space itineraries of tunneling orbits are related with the amounts of imaginary 
parts of actions gained by the orbits, so that the systematic search of significant tunneling orbits 
becomes possible. 
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I. INTRODUCTION 

Tunneling is one of the most typical and important 
phenomena in quantum physics, and for the past several 
years there is growing interest in natures of tunneling pro- 
cesses inherent in multi-dimensions. Quantum properties 
in multi-dimensional systems have been investigated ex- 
tensively in terms of classical dynamical concepts in the 
field of quantum chaos []], where the role of chaos, which 
is a generic property in multi-dimensional classical sys- 
tems, has been elucidated. It was found that tunneling 
properties are also strongly influenced by whether under- 
lying classical dynamics is chaotic or not 
, though tunneling process has no classical counterpart. 

Tunneling occurs typically between classical invariant 
components separated in phase space, such as between 
regular tori or chaotic seas. On one hand, mechanism of 
tunneling between distinct tori separated by chaotic seas 
has been studied in the context of chaos- assisted tunnel- 
ing 4], and its semiquantum analysis has been done, in 
which the diffusion process in the chaotic sea accompa- 
nied with tunneling paths from and into the tori around 
the boundaries of the sea is considered to dominate the 
tunneling transport [|| . Experiments have also been per- 
formed by measuring microwave spectra in the supercon- 
ducting cavity 0, and measuring momentum distribu- 
tions of cold atoms in an amplitude-modulated standing 
wave of light |l(J, |llj . 

On the other hand, tunneling between two chaotic seas 
separated by an energy barrier has been studied by sym- 
metric double wells p| ■ It has been shown that the spec- 
tra of tunnel splittings are reproduced by the orbits which 
consist of instanton process under the barrier and homo- 
clinic exploration in each chaotic well. 

Generic aspects of the link between tunneling pro- 
cess and real domain process in non-integrable systems 
have been examined in oscillatory scattering systems • 
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They made an energy domain analysis for a model with 
continuous flow, while in the present study we make a 
time domain one for a scattering map. The semiclassi- 
cal interpretation of complicated wave functions has been 
given in terms of oscillations of the stable manifold and 
an inherent property in flow systems, the divergent be- 
havior of movable singularities of classical solutions on 
the complex time plane. 

In near-integrable regime, the role of resonances has 
been elucidated in the tunneling transport between sym- 
metric tori, by means of classical and quantum pertur- 
bation theories [l3j. 

In any case, if one wants to know mechanism of tun- 
neling in chaotic systems by relating it with the under- 
lying classical structure, the use of complex orbits is in- 
evitable [T3 |. since tunneling is a purely quantum me- 
chanical process and is not describable in terms of real 
classical dynamics. Full account of such a process should, 
therefore, be given by complex classical dynamics. An 
attempt to make a full complex semiclassical analysis us- 
ing the complex classical dynamics has been performed 
to understand which kinds of complex trajectories de- 
scribe characteristic features of tunneling in the presence 
of chaos, and how the complex classical dynamics actu- 
ally enters into real physical process [5|, Il5l 1 1 6| . 

In Ref. || , it was found that the initial values of orbits 
which play a scmiclassically primary role form chain-like 
structures on an initial-value plane. A phenomenology 
describing tunneling in the presence of chaos based on 
such structures has been developed. 

In Ref. [15|, the first evidence has been reported 
which demonstrates the crucial role of complex-phase- 
space chaos in the description of tunneling process by 
analyzing a simple scattering map. Also in this case, it 
was found that initial values of orbits playing a semi- 
classically primary role form chain-like structures on the 
initial- value plane. 

Very recently, the chain-like structures are shown to 
be closel y re lated to the Julia set in complex dynamical 
systems The Julia set is defined as the boundary 

between the orbits which diverge to infinity and those 
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which are bound for an indefinite time. Chaos occurs 
only on the Julia set [l^. In Ref. 17], it was proven 
that a class of orbits which potentially contribute to a 
semiclassical wave function is identified as the Julia set. 
It was also shown that the transitivity of dynamics and 
high density of trajectories on the Julia set characterize 
chaotic tunneling. 

However, there still remains a problem in complex 
semiclassical descriptions. Significant tunneling orbits 
are always characterized by a property that the amount 
of imaginary parts of classical actions gained by the orbits 
are minimal among the whole candidates. It is, however, 
difficult to find such significant orbits out of the candi- 
dates, because an exponential increase of the number of 
candidates with time preventing us from evaluating the 
amount of imaginary part of action for every candidate. 

To solve this problem, in this paper, we investigate the 
structure of complex phase space for a scattering map 
model, and relate the structure to the amounts of imag- 
inary parts of actions gained by tunneling orbits. Our 
main idea is to relate the symbolic dynamics of a ho- 
moclinic tanglement emerging in complex domain to the 
behavior of tunneling orbits. It enables us to estimate 
the amounts of imaginary parts of actions gained by the 
tunneling orbits from symbolic sequences. 

The organization of the paper is as follows. In Sec. [HI 
the symbolic description of tunneling orbits is developed. 
This description requires an effective symbolic dynamics 
constructed on a complex homoclinic tanglement. In that 
section, we emphasize the importance of the application 
of the symbolic dynamics to tunneling, and the details of 
how we construct the symbolic dynamics itself is deferred 
to Sec. IIIII So it should be noted that in Sec. [H]we use 
the result in Sec. IIIII without any technical details. 

More precisely, in Sec. [n] the tunneling process is in- 
vestigated by the time-domain approach of complex semi- 
classical method. We introduce a scattering map which 
would be the simplest possible map modeling the energy 
barrier tunneling in more than one degree of freedom. 
Though real-domain chaos is absent in this model, it is 
shown that tunneling wave functions exhibit the features 
which have been observed in the systems creating real- 
domain chaos, such as the existence of plateaus and cliffs 
in the tunneling amplitudes and erratic oscillations on 
the plateaus. 

It is elucidated that such tunneling features originate 
from chaotic classical dynamics in the complex domain, 
in other words, the emergence of homoclinic tanglement 
in the complex domain. The symbolic description of the 
tanglement is introduced, and is applied to the symbolic 
encoding of the behaviors exhibited by semiclassical can- 
didate orbits. The amounts of imaginary parts of actions 
gained by the orbits are estimated in terms of symbolic 
sequences assigned to the orbits. Significant tunneling 
orbits are selected according to the estimation. 

Finally, tunneling wave functions are reproduced in 
terms of such significant orbits, and the characteristic 
features appearing in tunneling amplitudes are explained 



by the interference among such significant orbits. 

In Sec. IIIII the technical aspects which are skipped 
in Sec. |H] are described in full details. We first inves- 
tigate the construction of a partition of complex phase 
space. The homoclinic points are encoded into symbolic 
sequences by means of the partition. Then some numeri- 
cal observations are presented which relates the symbolic 
sequences and the locations of homoclinic points in phase 
space. On the basis of the observations, we study the re- 
lation between the symbolic sequences and the amounts 
of imaginary parts of actions gained by homoclinic or- 
bits. As a result, a symbolic formula for the estimation 
of imaginary parts of actions is derived. 

In Sec. IIVI we first conclude our present study, then 
discuss the role of complex-domain chaos played in semi- 
classical descriptions of tunneling in non-integrable sys- 
tems. It is suggested that the complex-domain chaos 
plays an important role in a wide range of tunneling phe- 
nomena of non-integrable systems. Finally some future 
problems are presented. 

II. SEMICLASSICAL STUDY ON TUNNELING 
PROCESS VIA COMPLEX-DOMAIN CHAOS 

A. Tunneling in a Simple Scattering System 

We introduce a simple scattering map model which will 
be used in our study. The Hamiltonian of our model is 
given as follows: 

+00 

H(q,p,t) = T(p) + V(q) Y, *(*-»). ( la ) 

Ti — — OO 

T{p) = p 2 /2, (lb) 
V{q) = fcexp(- 79 2 ), (lc) 

where k and 7 are some parameters with positive values, 
and the height and width of an energy barrier are given 
by k and 1 / ^f2rj respectively. A set of classical equations 
of motion is given as 

= f(<lj,Pj), ( 2a ) 

/:R 2 ^M 2 I (q,p) » (q + T'(p),p— V'(q + p)) (2b) 

where j denotes a time step with an integer value, and 
the prime denotes a differentiation with respect to the 
corresponding argument. K 2 denotes real phase space. 

Fig. shows the stable and unstable manifolds of a 
fixed point located at the origin (0,0). These manifolds 
are denoted as W s (0, 0) and W"(0,0) respectively. We 
should note that the map does not create chaos in real 
phase space, in contrast to the maps defined in bounded 
phase space such as the standard map |l9j| . One can rec- 
ognize this fact in several ways: for example, the map 
has only single periodic orbit, (q,p) — (0,0). It means 
that the topological entropy of the system is null. The 
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FIG. 1: Real-domain stable and unstable manifolds of the 
unstable fixed point located at the origin. 



other way is that W s (0,0) and W"(0,0) oscillate with- 
out creating homoclinic intersections. Any classical man- 
ifold initially put on the real phase space is stretched but 
not folded completely, and it goes away to infinity along 
W u (0,0). 

The quantum mechanical propagation for a single time 
step is given by the unitary operator: 



U = exp 



-V(q) 



exp 



-T{p) 



(3) 



where q and p denote the quantum operators correspond- 
ing to q and p respectively, which satisfy the uncertainty 
relation [q,p\ — ih. A quantum incident wave packet is 
put in an asymptotic region. The initial wave packet is 
given by a coherent state of the form: 



(q\q a ,p a ) 



i 



(7rfta 2 )V4 



exp 



(q-q a ) 2 .Pa{q a -^q) 



2ha 2 



21) 



(4) 

where a is a positive parameter and the width of the 
wave packet in the q direction is given by y/ha. q a and 
p a are the position and momentum of the center of mass, 
respectively. The initial kinetic energy is set to be far less 
than the potential barrier. 

Fig-IHshows the iteration of the quantum wavepacket. 
Several features are observed such as amplitude 
crossovers, the existence of plateaus and cliffs, and er- 
ratic oscillations on the plateaus. The same features have 
been reported in the case of dynamical tunneling in mixed 
phase space 5j . These are called the "plateau-cliff struc- 
ture", which has been confirmed in several systems as a 
typical structure of tunneling wave functions in the pres- 
ence of real-domain chaos |5J. However, as seen in our 
system, the existence of the plateau-cliff structure does 
not always need chaotic dynamics in real phase space. 
So, the features of wave functions observed here would 
be beyond our intuitive expectation based on the real 
classical dynamics. This strongly motivates the use of 
complex trajectories and complex semiclassical analysis 
to describe the features. 
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FIG. 2: (a) Squared amplitudes of wave functions, 
|(g|?7' l |g a ,p Q }| , evaluated quantum-mechanically for n — 2 
to 12 in every 2 time steps (ft = l, a = 10, k = 500, 7 = 0.005, 
q a — — 123, p a = 23). The incident wave packet is set in the 
asymptotic region in the side of q < 0. In the lower figures, 
the center of mass has been already reflected by the poten- 
tial barrier located around the origin, and those amplitudes 
observed in the transmitted region represent tunneling effect, 
(b) The enlarged figure at n = 10. 



B. Formulation of Semiclassical Analysis 

In the complex semiclassical analysis of our system, it 
is convenient to define a pair of canonical variables (Q, P) 
by 



Q 
p 



V2: 
a 

71 



-(p + iqa 2 ), 
(p- iq(T~ 2 ), 



and some notations by 

(Qo,Po) = {Q(qo,Po),P(qo,Po)) , 

(Qa,P a ) = {Q(q a ,Pa),P(q a ,p a )) , 



(5a) 
(5b) 

(6a) 
(6b) 



where (qo,Po) an d (qa,Pa) denote respectively the initial 
value of the classical map /, and the value specifying the 
center of the wave packet @- 

The wave function (q\U n \q a ,p a ) is represented by an 
n-fold multiple integral: 



(7) 



A n I dq .. .dq n—1 6Xp ~<S n . 
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which is a discrete analog of Feynman path integral, 
where 



Sn = Sn + Sq, 



(8a) 
(8b) 



S n = ]T [T(Pj-i) - V{qj)] (pj-i = qj - qj-i), (8c) 

3=1 

So = | [(Qo ~ Qa) 2 + {Pa + iQ a ){P a + iQ a - i2Q )] , 

(8d) 

A saddle point condition is imposed to the integral to give 
the semiclassical Van Vleck's formula, in which the wave- 
function can be expressed by purely classical-dynamical 
quantities. The saddle point evaluation with boundary 
conditions yields a shooting problem joining the initial 
and final points, (qo,po) and (q n -,p n )- In other words, 
the orbit should satisfy a set of classical equations of mo- 
tion and boundary conditions: 

(ij+uPj+i) = f (Qj,Pj) (0 < i < n), (9a) 
(qo,Po)eI, (9b) 
(OVuPn) e J 7 , (9c) 

where / : C 2 — > C 2 is the classical map extended into 
complex phase space, and X, T stand for manifolds de- 
fined by 

X = {(q, P )eC 2 I P(q,p) = P a }, (10a) 
T = {{q, P ) e C 2 I Jmq = 0}. (10b) 

We call the complex plane X the initial-value plane. 
The final condition in (|9c|l is required since we here want 
to see our wave function as a function of q n , which should 
take a real value. Therefore, the solutions of the classical 
equations of motion are given by a subset of the initial- 
value plane: 



M n = m/-"(f). 



(11) 



Since the initial "momentum" Pq is fixed as P a , the 
shooting problem will be solved by adjusting the initial 
"position" Qq in the initial-value plane X. 

The semiclassical Van Vleck's formula of the n-step 
wavefunction takes the form: 



{q n \U n \q a ,p a ) ~ A 



(qo,po)£M T , 



d 2 W n 



dq n dP 



exp - ( Sn 



The outline of the derivation of (|12|) follows the conven- 
tional one 20]. Further details are also given in the Ap- 
pendix of Ref. . 



C. Hierarchical Configuration of Initial Values 

As stated before, classical dynamics of the system 
does not create real-domain chaos. In contrast to that, 
the complex phase space has a complicated structure. 
Fig-EI a ) shows a typical pattern of Ai n , which consists of 
a huge number of strings. A single string covers the whole 
range (— oo, +00) of the final q n axis, in other words, each 
string has a solution of © for any choice of q n . So we 
call each string a branch. 

The morphology of M n is elusive. To describe it, we 
introduce the notions of chain-like structures and gener- 
ations. When blowing up any small area where branches 
are densely aggregated, one can always find the config- 
uration of branches as shown in Fig. I^b). The con- 
figuration is drawn schematically in Fig. |3fc), where a 
set of branches located in the center are linked to each 
other in the horizontal direction with narrow gaps, and 
form a chain-like structure [H Il5|. In any smaller area 
where branches are densely aggregated, one can find 
again chain-like structures with finer scale. 

As shown in Fig. Et c ) 5 small chain-like structures are 
arranged in both sides of the central large one, and 
the same arrangement repeats around each of the small 
chain- like structures. This observation means that the 
branches in M n have a hierarchical configuration. Then 
it may be natural to assign a notion of generation to 
each chain-like structure in the hierarchy. For example, 
in Fig. Etc), we can say that the first 4 generations are 
displayed. 

The hierarchical configuration of branches in M n can 
be understood more clearly by relating it to stable and 
unstable manifolds in complex domain. We will show two 
facts: first, chain-like structures are created by the itera- 
tion of a small area element along W s (0, 0) and W"(0, 0). 
Second, the intersection Ifl W s (0,0) forms the main 
frame of the hierarchical configuration of branches. In 
order to see these facts, we use the notation: 



(12) 



where the sum is over the complex orbits whose initial 
points are located on M. n just defined. <fi(qo,Po) is the 
Maslov index of each complex orbit. W n (q n , Pq) is a gen- 
erating function which yields a set of canonical transfor- 
mations as 



OWn 



dq„ 



= p n , and — - 

Po 9P 



(13) 



m = mr(o,o). 



(14) 



In the following, we mean a manifold with a real 1 di- 
mension (resp. real 2 dimensions) by a "curve" (resp. 
"surface" ) . 

When the classical map / is extended to the complex 
phase space C 2 , both W s (0, 0) and W n (0, 0) are surfaces 
in the complex domain at least locally. Since the initial- 
value plane X is also the same dimensional surface, the 
dimension of the intersection M is lower than 1 in gen- 
eral, i.e., the intersection is neither a set of curves nor 
a surface, but may be fractal like the Cantor set, whose 
Hausdorf dimension may be between and 1. Let (go, Po) 
be an element of M, and AX be a small neighbourhood 
of (qo,po) on the plane X. As the time proceeds, the orbit 
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FIG. 3: (a) M n on the initial-value plane X for n = 10. 
(b) The structure of A4 n found in any small area on X where 
branches are densely aggregated, (c) Schematic representa- 
tion of (b). Chain- like structures have a hierarchical config- 
uration in M n - A solid square located at the center of each 
chain- like structure represents an element of X n W s (0, 0). 



of (qo,Po) converges to the origin (0, 0) by definition. On 
the other hand, the iteration of AX approaches (0, 0) at 
the initial time stage, however, it is in turn spread along 
W"(0,0), and at a sufficiently large time step, it almost 
converges to W"(0,0). We note that this process is not 
specific to our system, but the one that the stable and 
unstable manifolds always have. Such a process takes 
place also in real-domain dynamics or even in integrable 
one. 




300 



(b) 













FIG. 4: (a) The n-step iteration of M n n AX for n = 10, pro- 
jected on real phase space. Bold and thin parts almost agree 
with W"(0, 0) in real and complex phase spaces respectively 
(in the thin part, those points which have quite large Imp„ 
are omitted), (b) Schematic representation of M„ n AX. The 
center dot represents a point in M. Hatched and non-hatched 
parts correspond to the bold and thin parts in (a) respectively. 




FIG. 5: M„ for (a) n = 2 and (b) n = 10, in the same 
range of the plane X. In (a), the central five branches forms 
a chain-like structure, which is the only chain-like structure 
at n = 2. Each filled square represents a point in M. As n 
increases, branches are densely aggregated in the neighbour- 
hood of each filled square as shown in (b). When enlarging 
the neighbourhood, one can always find the configuration of 
branches shown in Fig. ffl bl . 



According to this mechanism and the conditions in 
one may expect that for a sufficiently large time step n, 
the n-step iteration of A4 n C)AT almost agrees with a set: 



W"(0,0)n{( ? ,p)eC 2 | lmg = 0}. 



(15) 



Even in the case of n — 10, the n-step iteration of 
A4 n n AX, which is projected on the real phase space 
in Fig. Ufa), gives a sufficiently close manifold to that 
presented by (TSJ. 

The bold and thin parts represent the parts of the fi- 
nal manifold which almost coincide with the set (|15l) in 
real and complex phase spaces respectively. This means 
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FIG. 6: A variety of behaviors exhibited by orbits launching 
from M. Solid and broken lines represent Re q and Imq re- 
spectively. In (b) , the initial points of the trajectories belong 
to the 6th, 8th, and 10th generations respectively in the order 
from the top row to the bottom one. 



that the bold part have sufficiently small imaginary part 
of final momentum compared to that of the thin part. 
Fig. 0Jb) shows AX and the schematic representation of 
M n . The center dot represents (qo,po) appearing in the 
above process. Hatched part and non-hatched part corre- 
spond to the bold and thin parts in Fig.^Ja), respectively. 
Neighbouring branches are connected with each other via 
caustics defined by dq n /dQo\p = 0. The caustics are 
created by the oscillations of the real-domain W"(0, 0). 

In this way, the creation of a chain-like structure can 
be understood by considering the behavior of a small area 
clement first approaching real phase space with the guide 
of W s (0, 0), then spreading over W s (0, 0). In particular, 
the direction in which the real-domain W"(0, 0) stretches 
the area element determines the direction of the chain- 
like structure on the plane X. 

It was found that a chain-like structure on the plane 
X is always created around a point in M . Fig. [S] shows 
the creation of chain-like structures around the points in 
M as the time step proceeds. In Fig. Etc), the points in 
M are represented by filled squares. That is to say, the 
hierarchical configuration of chain-like structures implies 
those of the points in M . Accordingly, M constitutes the 
main frame of M. n , and the notion of generation is also 
assigned to each point in M. 

Since, as mentioned above, the iteration of initial 
points in chain-like structures to real phase space is de- 
scribed by orbits launching from M, the structure of such 
orbits will tell us that of orbits launching from A4 n , the 
latter structure is necessary for our semiclassical analysis. 
The study of orbits on the stable and unstable manifolds 
is suitable for more canonical arguments since they are 
compatible with the theory of dynamical systems |2l| . 

Fig.[S]shows a variety of itineraries of the orbits launch- 
ing from M. In Fig.^a), the top row shows a typical be- 
havior observed in M, where both Reg and Img oscillate 
in an erratic manner for some initial time steps and even- 
tually approach the origin. Regular itineraries such as 



periodic oscillations coexist among stochastic itineraries 
as shown in the middle row, where an approximately pe- 
riodic 2 behavior is seen. Another type of orbit is shown 
in the bottom row, where the trajectory first oscillates 
with period 2 and then turns into a periodic 3 motion. 
The close relation between itinerating behaviors and the 
notion of generation can be seen clearly in the case of 
periodic oscillations, as shown in Fig.|SJb). In each row 
of the figure, the period for which a trajectory oscillates 
agrees with the generation of the initial point of the tra- 
jectory. 



D. Emergence of Homoclinic Tanglement in 
Complex Phase Space 

Here we will show that the hierarchical structure of 
M, which is the main frame of Al n , is the manifestation 
of complex-domain chaos. First it is shown that homo- 
clinic tanglement emerges in the complex domain. The 
laminations of W s (0,0) and W tt (0, 0) are densely devel- 
oped around the origin (0, 0) in the complex domain, and 
they are described by the notion of generation. Then it 
is explained that the hierarchical structure of M is cre- 
ated as a consequence of the emergence of the complex 
homoclinic tanglement. 

In order to study the structure of complex phase space, 
we introduce coordinates on W s (0, 0) and W™(0, 0) as fol- 
lows. Let $ s and 4>„ be conjugation maps from a complex 
plane C to W s (0, 0) and W"(0, 0) respectively, which sat- 
isfy the relations: 

/.(&) = A^£ s for ^ e C, (16a) 
futtu) = \£u for^eC, (16b) 

where f s and f u are 

/.(£) = ($; 1 /*,)(0, (17a) 

and A denotes the maximal eigenvalue of the tangent map 
on the unstable fixed point (0,0) |22j. The above coor- 
dinates £ s and f„ are normalized in the sense that the 
classical map / is represented by a linear transformation 
on each coordinate. The set of homoclinic points asso- 
ciated with the origin, W s (0,0) n yV u (0,0), is obtained 
numerically on these coordinates as shown in Fig. [7| 

In each of Fig.[7fa) and (b), enlarging the neighbour- 
hood of the origin by A" times (n > 1), one obtains the 
same figure with the original one due to the relations in 
(|16ll . This means that the homoclinic points are accumu- 
lated around the origin (0,0), in other words, the lami- 
nations of yV u (0, 0) and W s (0, 0) in the complex domain 
are densely developed around the origin. It is a direct 
numerical evidence of the emergence of homoclinic tan- 
glement in complex phase space, and thus null topological 
entropy in real phase space does not always exclude the 
existence of chaos in complex domain. 
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FIG. 7: The set of homoclinic points associated with the 
unstable fixed point (0,0) plotted on (a) W s (0,0) and (b) 
W"(0, 0). The origin of each figure corresponds to the unsta- 
ble fixed point, and the horizontal axis through the origin is 
included in real phase space. 
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FIG. 8: (a) The intersection M plotted on the £ s plane, 
(b) The domain D plotted on the £ s plane with the set of 
homoclinic points superposed. 



Fig. Eta) shows M plotted on the £ s coordinate. Sim- 
ilarity to Fig. El^a) is evident, which suggests that the 
creation of the hierarchical configuration of M on the 
plane X is due to the emergence of homoclinic tangle- 
ment in the complex domain. The relation between the 
structures of M and the complex homoclinic tanglement 
is made more clear by the notion of generation, which we 
have already introduced. To see the relation, the precise 
definition of generation is given as follows. Let D be a 
connected domain in the £ s plane satisfying the condi- 
tions: 



(0,0) € D, 
fs{D) C D, 



(18a) 
(18b) 



where f s is the linear transformation defined in (|16af) . 
Then denoting D\f s (D) by D, the £ s plane is decom- 
posed into a family of disjoint domains as follows: 



U fs{D) 



C-{(0,0)}. 



(19a) 



(19b) 



Thus for any point t; s in this plane except (0,0), there 
exist a unique integer n such that £ s £ f~ n (D). We 



define the generation of the point £ s as such an integer 
n . 

The shape of the domain D is shown in Fig.[S{b). Note 
that the relations in 119fl hold irrespective of the shape of 
D. In our earlier publication 15] , D was chosen as a disk. 
In the present study, another choice of D is proposed. 
More precise definition of D is given in Sec. Ill El in the 
context of the construction of symbolic dynamics. 

For any point in the neighbourhood of the origin in D, 
the forward orbit approaches the origin straightforwardly 
at an exponential rate in the original coordinate (q,p). 
For any point in the n-th generation (n > 1), it takes at 
least n steps until the orbit starts to approach the origin 
exponentially, and thus it can exhibit a variety of be- 
haviors during its itinerary. That is why the oscillations 
shown in Fig. [5] are related with the generations. 

Fig. El shows the generations in real phase space. On 
one hand, in real-domain dynamics, a single iteration of 
the map / creates a single oscillation of W s (0, 0). On the 
other hand, by definition, any point of W s (0, 0) in any 
generation is mapped by / to another point in the neigh- 
bouring generation. Therefore in real phase space, each 
generation corresponds to a single oscillation of W S (Q, 0). 

In the complex domain, however, as seen in Fig.|3jc), 
higher generations describe finer scales of the hierarchical 
configuration of M on the plane X. Fig. ^c) also shows 
that for any point of M and in its small neighbourhood 
on the plane X, the number of the members in M which 
are included in individual generations increases exponen- 
tially with the generations. This situation clearly shows 
that the homoclinic tanglement of W s {0, 0) and W"(0, 0) 
creates the hierarchical configuration of M on the plane 
X. 

It should be noted that real-domain chaos, if it exits, 
also creates the hierarchical configuration. Fig. 1101 shows 
a simple analogy of our present situation with a horse 
shoe map on a two-dimensional plane. One-dimensional 
W s (0, 0) and W tI (0, 0) create homoclinic tanglement, and 
the one-dimensional bar across the tanglement corre- 
sponds to our initial-value plane X. In this case, both 
X n W s (0,0) and X n W u (0,0) form the Cantor sets, 
and the fractal structures of these intersections are orig- 
inated from densely developed laminations of W s (0, 0) 
and W u (0,0). The interval on W s (0,0) indicated by a 
bold curve corresponds to our domain D. Replacing C in 
(|19a|) by R, one can define generations in a similar way. 
Due to the horse shoe dynamics on this plane, W s (0,0) 
in the n-th generation has 2™ intersection points with the 
initial- value bar. 

In this sense, whether in real domain or in complex 
one, the hierarchical configuration of Ifl>V s (0, 0) on the 
plane X is nothing but the manifestation of chaos. In 
particular, the emergence of the above configuration in 
the complex domain is a piece of evidence for complex- 
domain chaos. Once we know that chaos exists also in 
the complex domain, the methodology studying chaos in 
the real domain can be applied to the analysis on the 
complex domain. In particular, symbolic dynamical de- 
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FIG. 9: Generations of W s (0,0) in real phase space. The 
bold and dotted curves represent even and odd generations 
respectively, except for the bold part emanating from the ori- 
gin, which represents the generations lower than or equal to 
zero. The generation increases monotonically as q — > — oo. 




FIG. 10: The horse-shoe map denned on a two-dimensional 
plane. The one-dimensional bar across the tanglement cor- 
responds to our initial-value plane T, and the bold curve on 
W(0, 0) corresponds to our domain D. The dotted line is the 
boundary of partition which creates binary codes. 



homoclinic orbits. We here present only the final re- 
sults of our study on this description, and the detailed 
explanation is given in Sec. IIIII On the basis of the re- 
sults, the symbolic description of semiclassical candidate 
orbits is discussed as follows. First the orbits launch- 
ing at M(=In W s (0, 0)) are encoded into symbolic se- 
quences. We compare both configurations of homoclinic 
points and the elements of M, on the £ s coordinate set on 
W s (0,0). A clear similarity between both configurations 
enables us to find, for each element of M, a homoclinic 
point located in the neighbourhood of the element of M. 
Since the homoclinic point is encoded into a symbolic 
sequence, we encode the element of M into this sym- 
bolic sequence. Then semiclassical candidate orbits are 
encoded into symbolic sequences. Since the behaviors of 
orbits launching at a single chain-like structure are de- 
scribed by the behavior of the orbit of an element of M 
located at the center of the chain-like structure, we as- 
sign the symbolic sequence of the element of M to the 
initial points of the chain-like structure. 

Symbolic dynamics is usually constructed by finding 
a generating partition Q, which is the partition of phase 
space satisfying the relation: 



-t-oc 



V no) 



(20) 



where the l.h.s. is the product of all partitions created 
by the iterations of Q, and the r.h.s. is the partition of 
phase space into its individual points |2l| (here we should 
consider the "phase space" as the closure of a set of ho- 
moclinic points). Roughly speaking, Q is the partition 
of phase space such that each separated component of 
phase space corresponds to a symbol and for every bi- 
infinitc sequence of symbols there may at most exist one 
trajectory of the original map. 

In our case, a partition of complex phase space is de- 
fined in terms of a phase part of the gradient of the po- 
tential function: 



scription of orbits, which is available if one finds a proper 
partition of phase space to define it, is a standard tech- 
nique in the theory of dynamical systems pH , and can be 
very useful tool to analyze complicated phase space struc- 
tures. Homoclinic orbits are also describable in terms 
of the symbolic dynamics, so our strategy to study the 
hierarchical configuration of M hereafter is to take the 
symbolic description of homoclinic orbits. 



E. Symbolic Description of Complex Orbits 

We will explain the symbolic description of the com- 
plex homoclinic orbits and its application to the symbolic 
description of semiclassical candidate orbits. The sym- 
bolic description of homoclinic orbits contains the con- 
struction of symbolic dynamics which works effectively 
and the estimation of imaginary parts of actions for the 



V'(q) = -2 1 kcxp[-A(q)-iB(q)], (21) 

where A(q), B(q) € R. The boundaries of our partition 
are defined as follows: 

{(?,P)£C 2 | B(q) = [2vxy-(3x+l)y/2]ir}, (22) 

where v is arbitrary positive integer, and (x, y) is an el- 
ement of the set: 

T = {(+!,+!), (-!,+!), (-1,-1), (+1,-1)}. (23) 

The intersection between the q plane and the boundaries 
of the partition are shown in Fig. ^] Each boundary is a 
three-dimensional manifold according to (|22|l . so that the 
intersection is a set of one-dimensional curves, x and y in 
(|22|l represent the signs of Re q and Im q of the points on 
a boundary respectively, i.e., the pair of x and y specifies 
the quadrant of the q plane where the q component of 



the boundary is included, v is the "distance" between 
the boundary and the origin (0, 0), in the sense that the 
axis Re q = Im q on the q plane intersects the boundary 
at q « (x-\/ utt/j, y^J 1/71-/7) f° r ^> 1- 

Then our partition, denoted by V, is defined as a set 
of phase space components as follows: 

V = {U(x,y,v) I {x,y,v)£S}, (24a) 
S = (T«)N)U {(0,0,0)}, (24b) 

where U(x, y, v) for (x, y, v) £ T®N, whose q component 
is displayed in Fig. ^2 denotes the region enclosed by 
two boundaries associated with (x, y, v) and (x, y, v + 1), 
and U(0, 0, 0) denotes the complement of the union of all 
components associated with (x, y, v) £ T®N. The origin 
of phase space, (0,0), is included in [7(0,0,0). 

Fig. 1121 shows the £ s plane divided by f n (P) for n = 
and 1. The central domain in Fig.^Ja) is the domain D, 
as is explained in Sec. Ill Dl The domain D is defined as a 
connected domain in U(0, 0, 0) n yV s (Q, 0) which includes 
the origin (0,0). If f n ('P) for any n > divides clearly 
the set of homoclinic points displayed here, then V is the 
generating partition, i.e., the relation l|20|) holds for V 
(it is sufficient to consider the case of n > since any 
homoclinic point is mapped to the region displayed in 
Fig. I12f a) by the iterations of the map). 

On one hand, / n ("P) for n = divides a set of homo- 
clinic points clearly as shown in Fig. I12f ah This means 
that our partition is a reasonable approximation of the 
generating partition Q. So by means of V 1 we can con- 
struct a symbolic dynamics which works effectively. On 
the other hand, there are some regions in the £ s plane 
where f n {V) for n > 1 fails to divide the set of homo- 
clinic points clearly. In such regions, it may be necessary 
to improve our partition to obtain the generating parti- 
tion, and it is our future problem. Note that our com- 
plex dynamics is not proven to be sufficiently chaotic, 
i.e., hyperbolic, and that the existence of the generating 
partition for non-hyperbolic systems is an open problem. 
The improvement of the partition is of a mathematical 
interest, however, as actually demonstrated below, the 
present definition of V is sufficient for our semiclassical 
analysis. 

In terms of the partition V, each homoclinic point is 
encoded into a bi-infinite symbolic sequence of the form: 

...OO a- n a_(„_!) . . . a_i . a Q a\ . . . a n O O . . . ,(25) 

where O = (0,0,0), n £ N, and a k £ S (\k\ < n). The 
symbol aft stands for a phase space component denoted 
by U(ak) which contains the fc-th iteration of the ho- 
moclinic point, whose initial location is in U(ao). The 
finite sequence of (|fc| < n) is accompanied with semi- 
infinite sequences of O's in both sides. It reflects that any 
homoclinic point approaches the origin of phase space by 
forward and backward iterations of the map /. In par- 
ticular, homoclinic points included in D are encoded into 
symbolic sequences of the form: 

. . .0 O a_„ a_(„_i) . . . a_i . O O . . . , (26) 
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FIG. 11: Intersection between the q plane and the 
boundaries of the partition. The curves labeled by a, b, c, 
and d indicate the boundaries associated with (x,y, v) = 
(1, 1, 1), (1, 1, 2), (1, 1, 3), and (1, -1, 1) respectively. 
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FIG. 12: Boundaries of (a) the partition V and (b) its forward 
iteration f(V) plotted on W s (0, 0), with homoclinic points su- 
perposed. The generations of the homoclinic points displayed 
are lower than or equal to 1. (c) The enlarged figure of the 
hatched part of (b). The boundaries of f(V) fails to divide 
the set of homoclinic points clearly, in the hatched part of 
(c). (d) The enlarged figure of the hatched part of (c). Ho- 
moclinic points are aggregated densely like a thick band. The 
boundary curves might touch the set of homoclinic points. 

where n £ N and afc £ S (0 < k < n). Any symbol in 
the r.h.s. of the decimal point is O, since the forward 
trajectory is always included in a component U(0) due 
to the definition of D. 

Further investigation of symbolic dynamics, which is 
presented in Sec. IIIII enables us to estimate the amount 
of imaginary parts of actions gained by homoclinic orbits. 
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Let w a homoclinic point which has a symbolic sequence 
of the form l|25|) . When we denote — {xk,Vk, v k) for 
k G Z, the imaginary part of action gained by the forward 
trajectory of w, denoted by s(w), is estimated as follows: 

+00 

s(w) = ^ImpTfo-O-VGfc)] 

~ ^ W 1/2 1/2 u I/ 2 1/2 \ 

(27) 

where 7 is a parameter of the potential function V(q). 

From these results, we discuss the symbolic description 
of semiclassical candidate orbits. To this end, we display 
schematically the configurations of homoclinic points and 
the elements of M, which are shown in the previous Figs. 
|7[ a) and IHfa) ■ The configurations can be seen clearly 
by introducing graphs on the £ s coordinate of W s (0,0), 
whose vertices represent the homoclinic points or the el- 
ements of M. 

The configuration of the homoclinic points in the do- 
main D is represented by a graph Koo defined as follows. 
First, by a graph K\ which is shown in Fig- ll^f a). we rep- 
resent homoclinic points which have symbolic sequences 
of the form: 

...O O a-x .0 O ... (a-i eS). (28) 

Second, we attach smaller copies of K\ to each vertex of 
the original K\ to obtain a graph K2 which is shown in 
Fig. I13f bl. By the graph K%, we represent homoclinic 
points which have symbolic sequences of the form: 

...O O a- 2 a-x . O O... (a_ 2 , a-i € S). (29) 

In a similar way, a graph K n is obtained for any integer 
n > 3 by attaching much smaller copies of K\ to each 
vertex of K n -x- Since inclusions hold as Kx C K2 C 
K3 C . . ., we finally obtain a graph which is shown 
in Fig. 1141 as the union of the graphs Kx, K2, K3, .... 
Koo is a tree graph, i.e., there is no loop on the graph. 

Next we consider the configuration of the elements of 
M. As shown in Figs, da) and Ha), there are clear 
similarities between configurations of homoclinic points 
and the elements of M. We define M n as a subset of M 
such that the members belong to generations lower than 
n. Then we observed numerically the following two facts: 
one is that on W s (0,0), the elements of M n for n > 3 
have the same configuration as a graph K n -i (M n for 
n = 2 has only a single element, and M n for n = 1 is 
null). The other is that on the plane X, the elements of 
M n are located at the centers of chain-like structures in 
M n (see Fig. He)). 

In the case of homoclinic points, each vertex of K n ^2 
corresponds to a symbolic sequence of the form: 

...OO a_(„_ 2 ) a_( n _ 3 ) ... a-x .O O .... (30) 





(a) 


\, (1.1.2) 




IjT, (1,1,1) \ 








mil S \i 




Ml* (1,-1,1) ^* 




ffl 1 U, -1, 2) 





FIG. 13: Graphs (a) K\ and (b) K2 representing ho- 
moclinic points which have symbolic sequences of forms 
...OOa-i . OO . . . and . . .O O a-2 a-i . OO . . . respec- 
tively, where a_2, a_i G S. In (a), the center vertex cor- 
responds to the case of a-i = O, and the vertices in the 
upper-right side and in the lower-right one correspond respec- 
tively to the cases of o_i = (1, 1, v) and o_i = (1, — 1, v) with 
v G N, and so on. 




FIG. 14: The graph Koa representing homoclinic points in- 
cluded in the domain D. Hierarchical configurations of ver- 
tices can be seen by enlarging the neighbourhood of any ver- 
tex. 

Then to each point of M n , one can assign formally a semi 
infinite symbolic sequence of the form: 

a-( n -2) a-{n-x) ■■■ a-x O O (31) 

Fig. 1151 shows some trajectories launching from M and 
symbolic sequences of the form (|31ll assigned to the ini- 
tial points of the trajectories. The signs and amplitudes 
of the q components at each time step are well described 
by (x,y)'s and z/s of the corresponding symbols. This 
means that there are clear similarities between motions 
of both trajectories launching from M and a set of homo- 
clinic points. More precisely, for an element of M which 
has a symbolic sequence of the form l|31f) . the forward 
trajectory is well approximated by that of a homoclinic 
point which has a symbolic sequence of the form: 

...OO. a_(„_ 2 ) a_(„_!) ... a-xO O .... (32) 
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FIG. 15: The orbits launching from M, and symbolic se- 
quences assigned to the initial points of the orbits. The up- 
per and lower signs in each symbol denotes x and y in (12 H 
respectively, and the integer denotes the value of v. A semi 
infinite part, O O O . . ., is omitted in each sequence. Solid 
and broken lines represent Reg and lmq respectively. The 
trajectories show (a) erratic motions, (b) approximately 2- 
periodic motions, and (c) monotonical approach to real phase 
space. 



In this way, the elements of M n are described by sym- 
bolic sequences. To the chain-like structure in M. n {<Z I) 
which is associated with each element of M n (c W s (0, 0)), 
we assign the same symbolic sequence as the element of 
M n , since as stated in Sec. Ill CI the motions of trajecto- 
ries launching from a single chain-like structure are well 
approximated, till they start to spread along the real do- 
main W"(0, 0), by the motion of a trajectory launching 
from an element of M located at the center of the chain- 
like structure. Then chain-like structures in M n are also 
described by symbolic sequences of the form 

Finally, we discuss the estimation of imaginary parts 
of actions gained by semiclassical candidate orbits. By 
making use of the above mentioned similarities between 
M and a set of homoclinic points, for an element of M„(c 
M) which has a symbolic sequence: 



a <zi a,2 . . . a„_ 3 O O . 



(33) 



(note that the notation is slightly different from (|31|I % 1 we 
evaluate the imaginary part of action of its orbit by a 
homoclinic orbit with a symbolic sequence: 



. O O . a,Q eti 0,2 



O O . 



(34) 



Due to this evaluation, the former imaginary part of ac- 
tion is estimated by applying i|27[l to the symbolic se- 
quence 

For each trajectory launching from a single chain-like 
structure in M. ni we approximate the imaginary part of 
action as that of the trajectory launching from the ele- 
ment of M n located at the center of the chain-like struc- 
ture. In this way, we can estimate the imaginary parts of 
actions gained by semiclassical candidate orbits. Though 
there are branches in A4 n which are not included in 



any chain-like structure, semiclassical contributions from 
them are negligible. See the discussion in Sec. Ill Fl 



F. Reproduction of Tunneling Wave Functions 

Semiclassical wave functions are constructed from 
significant tunneling orbits selected according to the 
amounts of imaginary parts of actions. The semiclas- 
sical mechanism of the tunneling process is explained by 
the structure of complex phase space. 

The wave function (q n \U n \q a ,p a ) given in @ is con- 
structed by two steps: First, out of M n (c M) the el- 
ements are picked up such that the imaginary parts of 
actions estimated by (|27[1 are minimal. Second, the Van 
Vleck's formula in (|12[) is evaluated for the initial points 
in chain-like structures associated with these elements 
of M n . Before picking up the elements out of M n , the 
Stokes phenomenon 23] should be taken account of. Fol- 
lowing the prescription given in Ref. ,2J| , we found that 
chain-like structures which have unphysical contributions 
to the wave function are associated with elements of M n 
whose symbolic sequences include symbols of the form 
or (—1,1,1/) with v € Z. At least, the chain- 
like structures which attain exponentially large semiclas- 
sical amplitudes are in this case. The justification is our 
future problem. See, the Discussion. 

After excluding such elements from M n , we have an 
ordering for elements of M n : 

w wi W2 ■■■ , (35) 

such that inequalities hold: 

< s(w ) < s{ Wl ) < s(w 2 ) < ■ . ■ , (36) 

where s(wk) for k — 0,1,2,... represents the imaginary 
part of action estimated by (|27|l for the forward trajec- 
tory of Wfc . 

From the Appendix II, it is seen that wq, which is 
primarily significant in M n for the wave function, has a 
symbolic sequence: 



ooo. 



(37) 



and the members of {wi, W2, ■ ■ ■ , w 2 ( n -2)}: which are sec- 
ondarily significant in M n , have symbolic sequences of 
the form: 



bbb...bbOOO. 



(38) 



where b = (1, 1, 1) or (-1, -1, 1). Due to J2TJ), the length 
of b b b . . . b b varies from 1 to n — 2, and the member 
of {wi , W2 , • • • , W2( n — 2) } which has b b b ... b b of length 
k — 2 (k < n) belongs to the fc-th generation. From l|27|) 
it is estimated that s(w\) — s(w2) = . . . = s(w2( n -2))- 

We observed that the orbits {wq, wi,...,W2( n -2)} 
behave as follows. The orbit u>o converges to the 
origin of phase space exponentially, so that it gains 
only small imaginary part of action. The orbits 
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FIG. 16: (a) The squared amplitudes of the quantum and 
semiclassical wave functions for n = 10, represented by dotted 
and solid curves respectively. Note that both curves are su- 
perposed, (b) Individual contributions from chain-like struc- 
tures to the semiclassical wave function shown in (a). For 
some components, symbolic sequences assigned to chain-like 
structures are presented, where b — (—1, —1, 1). 



{u>i, W2, ■ ■ ■ , W2(„-2)} first explore the vicinity of real 
phase space till the sequence b b b . . . b b terminates, and 
then converge to the origin exponentially. Such motions 
yield much smaller imaginary parts of actions than flip- 
ping motions in complex domain. The latters are generic 
trajectories launching from M . The motions of homo- 
clinic orbits whose symbolic sequences include sub se- 
quences of the form (|38|l are investigated in Sec. II I II 
There it is found that such homoclinic orbits explore the 
vicinity of real phase space. The motions of the orbits of 
wi, W2, ■ ■ ■ , W2( n -2) reflect those of homoclinic orbits. 

Fig. 1161 shows quantum and semiclassical wave func- 
tions for n = 10, the latter of which is constructed 
by taking account of the contributions from chain-like 
structures associated with wq, w% , . . . , u>2(n-2) • Both 
functions are in excellent agreement. The contributions 
from the other chain-like structures are much smaller 
than those from the chain-like structures associated with 
u>o, wx,..., u?2(n-2)- In fact, the squared amplitudes of 
them are of order ~ 10~ 50 at most. In particular, the 
contributions from trajectories which exhibit flipping or 
oscillatory motions are negligible, as shown in Fig. 1171 

On one hand, significant chain-like structures increase 
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FIG. 17: Squared amplitudes of the semiclassical components 
(the l.h.s.) which come from chain-like structures associated 
with the elements of M whose trajectories exhibit oscillatory 
motions (the r.h.s.). Solid and broken lines in the r.h.s. rep- 
resent Re q and Im q respectively. Flipping or oscillatory mo- 
tions in complex phase space gain large amounts of imaginary 
parts of actions due to large complex momenta. 



linearly with the time step n as seen above, and the 
branches included in individual chain-like structures also 
increase linearly with n, reflecting the oscillations of real 
domain W"(0, 0). So the numbers of branches which we 
should take into account increases algebraically with n. 
On the other hand, the number of chain-like structures in 
M n , and that of branches in M n increase exponentially 
with n. Therefore only a small number of branches are 
significant to describe the tunneling process. The alge- 
braic increase of the significant orbits is a consequence of 
the absence of real-domain chaos. When the real domain 
is chaotic, a small piece of manifold which has approached 
the real domain are stretched and folded without gaining 
additional imaginary part of action, so that the number 
of significant orbits can increase exponentially with time. 

The contributions from branches not included in any 
chain-like structures are also negligible, such as the 
branches in Fig. Ufa) except the central five ones. For 
orbits launching from such branches in M. n , there is no 
orbit launching from M which guides them to real phase 
space within n time steps. This means that these orbits 
have large imaginary parts of momenta at the time step 
n, so that they gain sufficiently large amounts of imag- 
inary parts of actions. If the imaginary parts of actions 
are positively large, the contributions from the orbits are 
small enough to be negligible, or if they are negatively 
large, the contributions from the orbits are unphysical 
due to the Stokes phenomenon, so that the orbits should 
be excluded from the whole candidates. 

The excellent agreement between both quantum and 
semiclassical calculations enables us to interpret semi- 
classically the features of tunneling wave functions. 
Fig. I16f b) shows a decomposition of the semiclassical 
wavefunction in Fig. llbT a - ) into components which come 
from individual chain-like structures. It can clearly be 
seen that the contributions from many chain-like struc- 
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tures reproduce the crossover of amplitudes in the re- 
flected region and the staircase seen in the transmitted 
region. We found that erratic oscillations on each com- 
ponent are due to the interferences between branches 
included in a single chain-like structure. The symbolic 
sequences of the form l)38|l can have various lengths of 
b b b . . .b b. The lengths of this part correspond to the 
generations of Wi, W2, ■ ■•, i£>2(n-2)- Thus both the stair- 
case and the crossover of amplitudes are created by the 
interferences between chain-like structures belonging to 
different generations. In this way, the complicated tun- 
neling amplitude is explained semiclassically by the cre- 
ations of chain-like structures on the initial-value plane, 
and by the exponential increase of the number of chain- 
like structures with time (though linear increase for sig- 
nificant ones) , which is due to the emergence of complex 
homoclinic tanglement. 

The semiclassical mechanism of the tunneling process 
in our system is summarized as follows. Stable and un- 
stable manifolds of a real-domain unstable orbit create 
tanglement in complex domain. The initial manifold rep- 
resenting a quantum initial state is located in the tangle- 
ment, so that the intersection points between the ini- 
tial manifold and the stable manifold form a hierarchical 
structure on the initial manifold. The orbits launching 
from the neighbourhood of each intersection point are 
guided to real phase space by the stable manifold, and 
then spread over the unstable manifold. The number 
of the orbits guided to real phase space increases expo- 
nentially with time (though significant ones increase al- 
gebraically), reflecting the hierarchical structure formed 
by the intersection points on the initial manifold. Then 
the interferences between these orbits create complicated 
patterns in the tunneling amplitude. 



III. SYMBOLIC DESCRIPTION OF COMPLEX 
HOMOCLINIC STRUCTURE 

A. Construction of Partition of Phase Space 

Here we construct a partition of phase space which en- 
codes homoclinic points of the origin (0, 0) into symbolic 
sequences and defines a symbolic dynamics which works 
effectively. The behaviors of homoclinic orbits and the 
estimation of imaginary parts of actions for the orbits are 
presented in the following Sees. IIII Bl and IIII CI respec- 
tively. 

There has been extensive studies on the construction 
of generating partition in real phase space |'2fil| even in 
non-hyperbolic regimes j2||. In such real-domain stud- 
ies, the boundaries of generating partition are roughly 
approximated by a set of folding points created by the 
single-step iterations of flat manifolds. However, the ex- 
tension of such working principle to complex phase space 
is not obvious. 

In order to know the locations of the boundaries of the 
generating partition, we consider the iterations of flat 



manifolds given by the form {(q,p) G C 2 \p = po} for a 
complex value po. Fig. 1181 shows the iterations of some 
small pieces of such a flat manifold, which exhibit a vari- 
ety of behaviors depending on the initial locations of the 
small pieces. Since the magnitude and phase of the gra- 
dient V'(q) are controlled respectively by the functions 
A(q) and B(q) which appear in i|21|) . the differences of 
the behaviors shown in the figure mainly come from the 
differences of the values of A(q). The contour curves of 
A(q) and B(q) are shown in Fig. IT51 

In order to know the locations of the boundaries of 
the partition far from the origin (0,0), we introduce a 
coordinate (u, v) in the q plane: 

u= [{Rcq) 2 - (Imq) 2 ] /2, (39a) 
v = Rcq -Iraq. (39b) 

On this coordinate, we obtain the estimations: 

A(q) = 7 [(Re(7) 2 -(Im(7) 2 ] - log ]?| 

= 2 7 u + 0(log|u|) (v : fixed, |u| -> 00), (40a) 

B(q) — 2jRe q ■ Im q — Arg q 

= 2 7 v + (log H) (it: fixed, |u| -> 00), (40b) 

or 2" ( v + 0(\u\^ 1 ) (v : fixed, \u\ -> 00), (40c) 

where 7 is a parameter of V(q). The above estimations 
imply that far from the origin (0, 0) the magnitude and 
phase of V'(q) are controlled by u and v respectively. 

On one hand, in a region of phase space with u 3> 1, 
which includes the real asymptotic region, V'(q) almost 
vanishes so that the behavior of a small piece of mani- 
fold there, typically observed as that of mi in Fig. 1181 
is of a free motion. On the other hand, in a region 
of phase space with u <C — 1, any small rectangle on 
the (it,i>) coordinate, {q S C | u G [uq,uq + Ait], v G 
[vq,vq + Aw]}, with Mo <C — 1 is mapped by V'(q) ap- 
proximately in an annulus of radii 27&exp(— 2 7 ito) and 
2jk exp [— 27(«o + Ait)] on the q plane. In Fig. O a 
small piece of manifold 7712 is set as Ai> =77/7, which ap- 
proximately corresponds to a single period of the phase 
B(q), so that its final manifold looks like an annulus 
on the q plane (one radius is much smaller than the 
other). Thus when we put the u component of mi as 
[vq, vq + n7r/ 7 ] for an arbitrary integer n, the final man- 
ifold of 77i2 covers almost the same range of the q plane 
n times, i.e., looks like a superposition of n annuli on 
the q plane. The simplest way to distinguish each annu- 
lus would be to locate a boundary of partition along the 
contour curve v = VQ+mr/j on the q plane for any integer 
n. From this observation, one can expect that far from 
the origin (0,0), the boundaries of the partition should 
take the form {(q,p) G C 2 | v = Vq + nir/~f} with Vq and n 
being some real value and arbitrary integer respectively. 

Around the origin (0, 0), as shown in Fig.[2Ul the loca- 
tions of the boundaries can be estimated by the single- 
step backward iterations of flat manifolds evaluated in 
real phase space. The dotted lines represent a set of fold- 
ing points of the folded manifolds, each point of which is 
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FIG. 18: Iterations of small pieces of a flat manifold given 
by {(q,p) eC 2 p = p } for p = 100.0 + i 50.0, with contour 
curves of the u component of the (u, v) coordinate superposed. 
One of the small pieces, mi, is hardly deformed by the iter- 
ations of the map, and remains to be at almost the same 
momentum as initially given. The other one, m,2, expands 
over a wide range of phase space (a hatched region) after a 
few step iterations. 




-40 Re q 40 

FIG. 19: Contour curves of the functions (a) A(q) and (b) 
B(q). In (a), each curve is given by A(q) = {1 + log(27)} /2 + 
0.4n for n £ Z, which approaches the axes Reg = ilmg in 
both directions of the curve. In (b), each curve is given by 
B(q) = nn/6 for n S Z with the branch — tv < Argg < n. 
Bold curves in the side of Re q > correspond to the case of 
B{q) = 0, and those in the other side corresponds to the cases 
of B(q) = — 7r(Img > 0) and 7r(Img < 0). The bold curves 
intersect the real axis at q — ±l/y / 2~7. 



FIG. 20: Single-step backward iterations of flat manifolds 
given by {{q,p) £ C 2 j p — po} ranging from po = —75.0 to 
75.0 by every 25.0. The dotted lines represent a set of points 
each of which gives the maximal or minimal momentum on 
each folded manifold. 

defined by the condition dp(q)/dq — where the differ- 
entiation is along each folded manifold. The set of such 
folding points are given by the single-step backward it- 
erations of the lines {(q,p) € R 2 | q = ±1/^/2^}, so that 
the boundaries of the partition in complex domain are 
expected to intersect the real domain around these lines. 

We propose a partition of phase space which satisfies 
the rough estimations presented above both in complex 
and real domains. To define the partition, we prepare 
the notations: 

U(x,y) = {(q,p) G C 2 \xReq > 0, ylmq > 0}{41a) 
P{x,y,u) = [2uxy-(3x + l)y/2]w, (41b) 

where (x, y, v) is an element of T ® N which appeared in 
(I24bll . U(x,y) covers a single quadrant of the q plane and 
/3(x, y, v) always takes an integer times tt. For [x, y, v) G 
T ® N, a phase space component U(x, y, v) is defined by 

U(x,y,v) = {{q,p) | (q,p) G U(x,y), 

[B(q) - 0(x, y, i/)] [B(q) - 0(x, y,v+ 1)] < 0}, 

(42) 

and for (x, y, v) = (0, 0, 0), by 

U(x,y,v) = C 2 - (J U(x',y',p'). (43) 

(x' ,y' ,v')eT ®N 

Then our partition V is denned as a set of the above 
phase space components, as already shown in l|24a|) . and 
the boundaries of V take the form in Q22)). Such defi- 
nition of partition satisfies our rough estimation for the 
locations of the boundaries. In fact, in the complex do- 
main far from the origin (0,0), due to (|4U|) . the relation 
B(q) — /3(x,y,v) leads to v i=s vq + nixj^j when we set 
vq = — (3x + l)j/7r/47 and n — vxy. Moreover, Fig. lTyT b) 
shows that the boundaries of V indicated by B(q) = 
0(x, y, v) for (x, y, v) = (1, 1, 1), (1, -1, 1), (-1, 1, 1, ) and 
(—1,-1,1) intersect the real phase space at q = ±1/ y/2j. 



15 



B. Properties of Homoclinic Orbits 

By the partition constructed above, each homoclinic 
points is encoded into a sequence of symbols of the form 
(x, y, v). In order to estimate imaginary parts of actions 
gained by homoclinic orbits, it is necessary to under- 
stand typical behaviors exhibited by the homoclinic or- 
bits. Here we present such typical behaviors as two ob- 
servations obtained from numerical computations. The 
first observation is concerned with the relation between 
the integer v, which is a member of a symbol (x, y, v) in 
a symbolic sequence, and the flipping amplitude of the 
corresponding trajectory. The other of them is concerned 
with the relation between the length of a consecutive part 
b b ... b (b £ S) in a symbolic sequence and the behavior 
of the corresponding trajectory. These are numerical ob- 
servations and we have no mathematical proof, but the 
phase-space itinerary of any homoclinic orbit can be well 
explained by the combinations of the behaviors presented 
in these observations. 

Before presenting the observations, we see the single- 
step folding process of flat manifolds to estimate the 
locations of homoclinic points in phase space. Let rrii 
and m/ be two-dimensional flat manifolds defined by the 
conditions p = pi and p = pf respectively, with pi and 
Pf being complex values. The intersection f{rtii) Dirif 
is given by {(q,p f ) £ C 2 | A(q) = -\og\c\, B(q) = 
— argc ± 2vtt, v = 0, 1, 2, . . .} where A(q) and B(q) are 
the functions given in l|4TJI) . and c = (pf — pi)/(2jk). 
Since the q components of the intersection points on m/ 
are located on a single contour curve of A(q), they are 
located along the axes Re q = ±Im q asymptotically as 
v — > +oo. Also the q components of homoclinic points 
are located along these axes asymptotically as v — > +oo, 
where v is a member of a symbol (x,y, v) in a symbolic 
sequence. 

Observation 1. Let {wi, W2, W3, . . .} be a set of homo- 
clinic points whose symbolic sequences take the forms: 

. . . a_2 a_i .61 ai a 2 . . . for w\, 
. . . a_ 2 a_i . b 2 ai a 2 ■ ■ ■ for w 2 , 
. . . a_ 2 a_i . b 3 ai a 2 ■ ■ ■ for w 3 , 

(44) 

where (fc 7^ 0) is any member of iS defined in l|24|l . and 
b v = (x,y,v) for v = 1,2,3,... with (x,y)(£ T) being 
fixed. Then the following relations hold: 
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FIG. 21: The trajectories of homoclinic points w v 
for v — 1,2,..., 5. The left and right columns dis- 
play the real and imaginary parts of the trajectories 
respectively. The symbolic sequence of w v is given 
by . . . O O (-1, 1, 1) . K (1,1,2)000... with K being (a) 
(1, 1, v), (b) (-1, 1, v), (c) (-1, -1, v), and (d) (1, -1, v). The 
dotted lines represent the case that {>„ = O. 



The r.h.s. of each equation denotes a pair of signs of 
real and imaginary parts. 

Fig.l2"Tlshows the trajectories of w„'s for v = 1,2, ... ,5. 
In these figures, one can see two facts: one is that the 
signs of Rego(w i y) and Im qo(w„) are described respec- 
tively by x and y in the symbol b v . The other is that the 
amplitudes of Keqa(w u ) and Imgo(Wi,) increase with v 
much faster than the amplitudes at the other time steps. 
Due to the second fact, the following approximations hold 
for large z/s: 



Pa 
P-i 



q-i 



-qo, 
qo. 



(46a) 
(46b) 



lim v / 7/(* /7r ) qa{Wv) = 0>J/)> 



lim Vt/O^t) p Q {w u ) = -(x,y), 



lim VVC^tt) P-x{w v ) = (x,y), 



(45a) 
(45b) 
(45c) 



where (g , po) is the current location of w v in phase 
space, and p_i is the momentum at the last time step. 



so that the sign of Hepo(w u ) (resp. Impoiw^)) is oppo- 
site to that of Kep-i(w u ) (resp. lmp-i(w v )) for large 
i/s. Fig.l22lshows an(w,,). po(w u ), andp_i(u>„) for much 
larger z/s. The magnitudes of the real and imaginary 
parts of these quantities are shown to have the depen- 
dence of the form y/vir/j for sufficiently large ^'s. Fig. 1231 
shows that qo{wj) diverges much faster than q_i(wj) and 
q-2(wj). Similarly , qo(wj) diverges much faster than 
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Qk{wj) for any other k 0, though that is not displayed 
here. 

The relation (|45a|l leads to 

lim \Reqo(w v )/Imqo(wv)\=l. (47) 

v — >+oo 

In the following, we explain why the relations (|45|l should 
follow, assuming that the relation 147|) holds for the ho- 




v 



FIG. 22: The dependences of (a) qo(w v ), (b) po(w v ), and 
(c) p-i(wv) on the subscript number v. The symbolic se- 
quences of w v are given in (a), where A v , B v , and C v denote 
(1, 1, v), (— 1, 1, u), and (—1,-1,!/) respectively. The dotted 
line in each figure represents \J vk /-y. In (a), all curves al- 
most coincide. 

moclinic points w v (v = 1,2,3,...) given by l|44l) . and 
that qo(w u ) diverges much faster than qk(w v ) for any 
k 7^ as v — ► +00. 

The relation (|45a|l is explained as follows. Since 
b v = (x,y,f), w v is included in a phase space compo- 
nent U{x,y,v) defined in (@2J and ig^J. Then the v 
component of qo(w„) in the (u,v) coordinate diverges 
as v — ► +00, since B(qo(w v )) diverges as v — > +00 
due to (22), and v « B(<fo(uv))/(27) « 2uxyn/(2j) 
for large i/s due to (gUJ) and Therefore from (|4"7|) 

and the relation Re go 1 Imtfti = v ~ vxyn/j, we obtain 
go(wv) ~ y/ ^71 /7 (a; y) for large z/s. 

The relation l|45bj) is explained as follows. The classical 
equations of motions in lead to the relation: 

qo(wu) + q- 2 {w v ) = 2q_ 1 (w v )-V'{q-i{w v )).{4&) 

Since qo{w l/ ) diverges much faster than q_2(w„) as v — > 
+00 due to our assumption, the r.h.s. of the above re- 
lation diverges in this limit. It means that q-\{w u ) also 
diverges as v — > +00, since the r.h.s. of the relation is an 
entire function of q-\. In particular, the u component of 
q-i(w lJ ) diverges as v — ► +00, since (q—i(w u ),p-i(w u )) 
is always included in the phase-space component U(a-i) 
irrespective of v. Furthermore, the u component di- 
verges to —00, since if the u component diverges to +00 
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that as the length of b b . . .b increases, the location 
of the turning point diverges, so that the trajectory 
corresponding to b b ... b does not approach to any point 
in phase space in the limit of the length. 

Observation 2. Let { wi, W2, W3 . . .} be a set of homoclinic 
points such that the symbolic sequences take the forms: 

. . . a_2 a_i . 6 6 ao ai 02 . . • for w±, 

. . . a_2 ei-i . b b b b ao <xi aa • • • f° r W2, 

. . . a_2 a_i . b b b b b b ao ai a% . . . for u>3, 

; (49) 

where b ^ (0, 0, 0) and the length of b b ... b for Wj (j — 
1, 2,3,.. .) is 2j. Then the trajectory of wj correspond- 
ing to the consecutive part b b ... b is included in a single 
phase-space component U(b), and the momentum almost 
vanishes at time step j — 1. Moreover, the following in- 
equalities hold: 

< Reqk-i(wj)/Reqk[wj) < 1 (0 < fc < j), (50a) 

< lm.q k {w j )/~hnq k - 1 (wj) < 1 (0 < k < j), (50b) 

< Req k (wj)/Reqk-i(wj) < 1 {j < k < 2j), (50c) 

< Imq k -i(wj)/Imq k (wj) < 1 (j < k < 2j). (50d) 



FIG. 23: The v dependence of (a) |<3 , -i(w^)/(?o(w I y)| and (b) 
\q~2{w v ) / q-i(w v )\ (the vertical axes of (b) and the inset have 
the common logarithmic scale) . The symbolic sequences of w v 
are the same as in Fig. 1221 Solid lines are the guides for the 
eye. 

with (q_i(uv),p_i(uv)) being in the same phase-space 
component, then g„i(w„) approaches the real axis and 
V'(q-i(wv)) vanishes, so that Im [q^Wv) + q-2{w v )] — 
Im [2<7_i(w„) — V (q-i(w u ))\ — > 0. However, the above 
contradicts that as v — > +00, |Im[(7o(w„) + q-2{w v )\ \ ~ 
llmgo^i/)! ~ \J VTfpi —* +00. Thus the u com- 
ponent of q-\{w u ) diverges to —00 as v — > +00. 
When the u component of q_i(w v ) is negatively large, 
V (M-i( w v)) is exponentially larger than q_i(w ll ), so 
that qoiwv) « —V'(q-i{w v )). This relation means that 
I g-i I ~ \/7 _1 log \qo\, and thus Pq{w v ) has the same de- 
pendence as qo{w„) on v due to the relation po = go - 
The relation l|45c|) is explained in a similar way. 

We proceed to the other observation. In usual sym- 
bolic dynamics, a consecutive part bb...b of a single 
symbol b in a symbolic sequence always corresponds to 
a fixed point in phase space or the motion approaching 
the fixed point. However, our classical dynamics always 
has only a single fixed point (0, 0) for any choice of 
positive parameters k and 7, which is easily checked 
by solving f(q,p) — (q,p), so that the phase-space 
motion corresponding to a consecutive part b b . . .b with 
b ^ (0,0,0) is not obvious. Our second observation says 
that the phase-space motion corresponding to the above 
consecutive part has a turning point. It is conjectured 



This observation is exemplified in Fig. El In the case 
that the length of b b . . .b in H49J) is given by 2j + 1 for 
Wj , (|50a|l and i|50b|l hold in the range of < k < j, and 
H50c|) and i|50d(l hold in the range of j < k < 2j. In this 
case, the momentum p k (wj) at k = j is quite small, but 
does not vanish. 

We conjecture that the q component of the turning 
point, qj(wj), diverges with the length of bb...b, i.e., 
the following relation holds: 

lim qj{wj) = (x-oo, 0), (51) 

where x is the sign of the infinity, which is given by the 
member of the symbol b = (x,y,u). This conjecture is 
based on the following observation. 

Fig. [23 shows the locations of the turning points qj{wj) 
ranging from j = 1 to 100, and the magnitudes and 
phases of the potential function V(q) at the turning 
points. It can be seen that as j increases, the magni- 
tude of V(qj(wjj) decreases to zero algebraically, and 
the phase of V(qj(wj)) approaches some constant. That 
is, V(qj(wj)) sa aj~^e ie for large j's, where a, (3 (both 
positive), and 8 depend on a symbol b which consists of 
a consecutive part b b ... b in a symbolic sequence of Wj . 
If this approximation holds for any large j, solving the 
equation V(qj(tVj)) = aj~@e lS , one obtains the solution: 

Uj ( Wj ) = D81ogi-log(a/*)]/(27), (52a) 
Vj ( Wj ) = -0/(2 7 ), (52b) 
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FIG. 25: (a) The locations of the turning points, qj(wj), 
for the trajectories of Wj's ranging from j = 1 to 100. The 
symbolic sequences of Wj's are the same as is given in the 
caption of Fig. 1241 (b) The magnitude and phase of the 
potential function V(q) at the turning point qj(wj) for ranging 
from j = 1 to 100. The phase, — 2 r yYieqj ■ Im§j, is plotted 
without taking mod 2n. In (a) and (b), solid lines are the 
guides for the eye. 



FIG. 24: (a), (b) Trajectories of Wj's for j = 5, 10, and 
15. The real and imaginary parts of the trajectories are dis- 
played in (a) and (b) respectively. The symbolic sequence of 
Wj is given by... O O . b b . . . b O O where the symbol 
b is given in (b) and the length of bb. ..b is 2j. The real 
and imaginary parts of each trajectory have the maximal and 
minimal amplitudes respectively, at time steps j — 1 and j. 
That is, the momentum almost vanishes at time step j — (c) 
Trajectories of Wj's with the same symbolic sequences as in 
(a) and (b) but for j = 100 and b = (1, 1, v) with v = 1, 2, 3. 
The vector field of V'(q) is superposed with the length nor- 
malized. In (a),(b), and (c), solid lines connecting points of 
trajectories are the guides for the eye. 

where (uj(vjj),Vj(wj)) is the location of qj{wj) on the 
(u, v) coordinate defined in This solution suggests 

that Heqj(wj) diverges and Im. qj(wj) vanishes as j — > 
+oo. Moreover, according to the Observation 2, qj(wj) 
and the q component of U(b) are included in the same 
quadrant of the q plane. Therefore the relation l|51|l is 



obtained. The justification of this relation needs further 
investigation of classical dynamics, and we hope to report 
the result of this issue elsewhere. 

We have shown that there are two types of behaviors 
exhibited by homoclinic orbits. In our numerical com- 
putations, the behavior of any homoclinic orbits can be 
understood by the combinations of only two types of mo- 
tions, one of which is the flipping motions almost along 
the axes Reg = ±Imq, and the other of which is the mo- 
tions almost along the contour curves of the v component 
in the (u, v) coordinate. Which type of motion occurs 
in the process from (qk,Pk) to (q k +i,Pk+i) along a single 
homoclinic trajectory depends on whether the neighbour- 
ing symbols in a symbolic sequence, a k and dfc+i, are the 
same (the latter type) or not (the former type) . The for- 
mer type of motion is characterized by the Observation 
1, and the latter one by the Observation 2. 
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C. Estimation of Imaginary Parts of Actions 

The imaginary parts of actions gained by the orbits 
of homoclinic points are estimated from the symbolic se- 
quences assigned to the points. We first consider the 
homoclinic points appearing in the Observations 1 and 
2, and the estimations of imaginary parts of actions for 
these cases are presented as the Observations 3 and 4 
respectively. Then using the latter two observations, we 
estimate the imaginary part of action for any homoclinic 
trajectory. The Observation 3 says that the imaginary 
part of action diverges linearly as v — > +oo, where v is 
a member of a symbol (x, y, v) in a symbolic sequence. 
The Observation 4 says that the amount of the imagi- 
nary part of action is bounded even if the length of a 
consecutive part b b ... b (b S S) in a symbolic sequence 
tends to infinity. In particular, we observed that the 
itinerary described by b b ... b gains little imaginary part 
of action compared to the other itineraries in a trajec- 
tory. This means that the homoclinic orbits appearing 
in the Observation 4 can play a semiclassically significant 
role. The Observations 3 and 4 are also entirely based 
on our numerical computations, and so far, we have no 
mathematical proof for these observations. 

For the trajectory of any homoclinic point w, here we 
consider the following imaginary part of action: 

s(w) = £)lm[T(p fc _i(ti;)) - V(q k (w))} , (53) 
fc=i 

where (qk(w),pk{w)) denotes the k step iteration of 
w. The sum in the r.h.s. converges due to the ex- 
ponential convergence of the trajectory to the origin 
(0,0), and it is the long-time limit of ImSVi which 
is given in (|8c|l . In the definition of s(w), we only 
take account of the contributions from the forward 
trajectories, since semiclassical wave functions in the 
time domain are determined by them. We do not 
define s(w) as the long-time limit of ImSVi which is 
given in JSEJ, since it includes an additional term So 
which depends only on the choice of an incident wave 
packet, not on classical dynamics. First we estimate s(w) 
for the homoclinic points appearing in the Observation 1. 

Observation 3. Let {wi,W2,wa, . . .} be a set of homo- 
clinic points which appears in the Observation 1. Then 
for any integer n > 1, the following relation holds: 

km [7/(2w)]i(/-»K)) = xy. (54) 
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FIG. 26: (a) The absolute values of the imaginary parts 
of actions, \s(f~ n (w v ))\, for n = 4 and v = 1 to 1000. 
The symbolic sequences of w v 's axe the same as in Fig. 
1221 The dotted line represents 2vkI^. (b) The abso- 
lute values of the kinetic parts and the potential ones of 
s(f~ n (w v )). The upper dotted line and the lower dotted 
curve represent 2vn and 2/(^/2^71-/7) respectively, where 
f{x) = :r/(2v / 7lcg^). (c) |ImT(p-i(xo„)) + ImT(p (w v ))\, 
\ImViq-i(w v ))+ImViqi(w v ))l and \hn.V(q (Wu))\. The top 
dotted line and middle dotted curve are the same as in (b), 
and the bottom dotted curve represents { 27 2 /( -y/22/ 71-/7) 
In each figure, the solid lines are guides for the eye. 



Fig-HU a ) shows the magnitudes of s(/ _n (iy_,)) for n = 4 
and v = 1 to 1000. It can be seen that \s(f~ n (w v ))\ w 
2vTi/ r y for large v's. The condition n > 1 in l|54(l is im- 
posed to always take account of the contributions from 
the flipping motions from q_i(uv) to qo(w u ) displayed in 
Fig. ED 



In the following, we explain the relation 1)54(1 . assum- 
ing that the Observation 1 holds, and that for the ho- 
moclinic points appearing in the Observation 1, qk(w u ) 
(resp. g_fc(w y )) for k = 0, 1, 2, . . . diverges much faster 
than q k+ i(w v ) (resp. q_( k+1 )(w v )) as v -> +00. Fig. [23 
suggests that the second assumption is valid for k — 
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and 1. For the other fe's, it has not been found numeri- 
cally whether the assumption is valid or not, since q k {w v ) 
and q-k{w u ) for k > 2 remain to be immediate values 
even for v k, 1000, so that the numerical computation 
needs too high accuracy to make clear the asymptotic 
behaviors of \q k \ and \q~k\ with sufficiently large magni- 
tudes. However, the exponential dependence of V'(q) on 
u and v shown in (14011 means that the large difference in 
1k(wj) (resp. q_ k )(wj)) results from the slight difference 
in qk+iiwj) (resp. q_ {k+1) (wj)) by the map f- 1 (resp. 
/), so that the assumption for k > 2 is expected to be 
valid. 

The relation (|54l) is explained as follows. In the kinetic 
part of s(f- n (w v )), i.e., E+™_ n+1 ImT(p k -i(w v )), each 
term is written as 

lm.T(p k ^i{w l ,)) = \Repk-\{w v )] \Lmp k -i{w v )] . (55) 

Due to the assumptions we put, the following inequalities 
hold for k ^ 0, — 1 and for large z/'s: 

\^Bp k {w v )\ < |Rep_i(w„)|, |Rep (wv)|, (56a) 
\lmp k (w v )\ < \Imp^i(w l/ )\,\Imp (w u )\. (56b) 

Then the kinetic part of s(f~ n (wi / )) is dominated by 
the terms Im T(p_i(-uv)) and lm.T(p {w v )) due to 
Since the Observation 1 says that the quantities in the 
r.h.s. of the above inequalities have v 1 / 2 dependences 
on large f's, ImT(p_i(w„)) and ImT(po(wv)) have lin- 
ear dependences on large v's due to (|55|l . Therefore the 
kinetic part of s{f~ n (w v )) is expected to has a linear de- 
pendence on large z/'s. Fig.l26lshows that the kinetic part 
of s(/~"(w„)) is actually dominated by ImT(p_i(u;„)) 
and 1thT{pq{w u )), and has a linear dependence on large 
z/'s. 

In the potential part of s(/~ n (w t/ )), i.e., 
T,^^_ n+1 lm.V(q k (w v )), each term is written as 



lu\V{q k {w v )) 



Im 



q k +i(w v ) - 2q k (w v ) + q k -i(w v ) 



2~/q k (w„) 



(57) 



by incorporating the classical equations of motions given 
in @ with the relation V'(q) = —2^qV(q) satisfied by 
our potential function. A simple calculation of the r.h.s. 
of ltST|) yields an inequality: 



\ImV(q k (w v ))\ < 



\q k +i(w v )\ + |gfc_i(w„)| 
2j\q k (w u )\ 



(58) 



Based on the assumptions imposed here, we can 
present the same discussion as that below l|47[) (note 
that the assumptions here are stronger than those im- 
posed there). As a result, one obtains the relations 
\qi(w^)\,\q-i(w u )\ w y /-f- 1 \og\q {w„)\ for large v. In 
a similar way, one obtains the relations \q±(k+i)( w u)\ ~ 

VT IT logT&fe(«V)l for l ar g e v ' s and k = 1,2,3,.... By 



using these relations, the r.h.s. of l|58[l is approximated 
by 

[27 2 /(k fe K)l)] _1 

f{\lk-i{w v )\) 
f(\q k+1 (w u )\) 



for 


k = 0, 


(59a) 


for 


k > 0, 


(59b) 


for 


k < 0, 


(59c) 



where fix) 



x/(2y/j log a;). Here we used an approxi- 
\q\ for large 



mation \q\ + y/^y 1 log |g| 

Since 1 90(^)1 is approximated as yj2vTi/ r ) for large 
z/'s according to the Observation 1, the r.h.s. of l|58() is 
approximated by 



2 7 2 / ^ 



/ (/2W7 
/(^/log'(log'(...log'(log'(27rz,/ 7 ))...)) 



for jfe = 0, (60a) 
for jfe = ±lj(60b) 

for \k\ > 2, (60c) 



where log's = (27)^ log ir, and the argument of the 
square root in (|60c|l is a \k\ — 1 fold logarithm of 2wv/j. 

Since fix) in (|60|l is monotonically increasing for large 
x, one can expect that the potential part of s(f~ n (w v )) 
for large z/'s is dominated by the terms ImV(q-i(w v )) 
and lmV(qi(w„)). More precisely, from (|60b|l . the po- 
tential part of s(f n (w u )) is expected to be approxi- 
mated by 2f{yj2vir /j) for large z/'s. And also, from 
(|60a|l . Im V(qo(w„)) is expected to vanish as v — > +00. 
Fig. shows that the potential part of s(f^ n (w l/ )) 
is actually dominated by the terms Im7(g_i(w„)) and 
Im V(qi(w v )) for large z/'s, and the asymptotic behav- 
ior o f the p otential part for large z/'s is described by 
2f(\j2vir /j). It is also shown that ImV (qo{w v )) tends 
to vanish as v increases. 

Consequently, due to the v dependencies of kinetic and 
potential parts as seen above, s(/~"(w 1/ )) for large z/'s is 
dominated by the kinetic part, which has a linear depen- 
dence of large z/'s. Therefore the following estimation is 
obtained for large z/'s: 

s(/-"K)) « ImT(p_iK))+ImT(poK)) 
= Rep_i(uv) • Imp_i (w u ) 

+Rep (w u ) ■ Impo(^) 
« (2i'TT/ r y)xy. (61) 

In the last approximation, the relations in l|45|) are used. 

Next we estimate s(w) for the homoclinic points 
appearing in the Observation 2. 

Observation 4. Let {wi,W2,w%, . . .} be a set of homo- 
clinic points which appears in the Observation 2. Then 
for any integer n, a sequence of imaginary parts of ac- 
tions, 



s(f- n (w 1 )) s(f- n (w 2 )) s(r n (w 3 )). 
is bounded. 



(62) 
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Fig. Eta) shows s(f~ n (w 3 )) for n = 4 and j = 1 to 100. 
It can be seen that s(f~ n (wj)) (j = 1, 2, 3, . . .) does not 
deviate largely from s(f~ n (w)) for a homoclinic point w 
which has a symbolic sequence: 



OO.b OO. 



(63) 



In the case of Fig. Eta), s(/~"K)) (j = 1,2,3,...) 
is mainly gained by the flippings of trajectories between 
phase-space components U(b) and U(0) (O — (0,0,0)). 
Though not displayed here, also in the case of the sym- 
bolic sequence: 



. a_2 o,—i .6 b . . .b ao a\ 



(64) 



with . . . a_2 a_i 7^ . . . O O or ao ax ■ • ■ 7^ OO..., the 
imaginary part of action does not deviate largely from 
that in the case of the sequence: 



. a_2 d-i -b ao ai 



(65) 



Therefore the imaginary parts of actions are mainly 
gained by phase-space itineraries described by sequences 
other than b b . . .b. This means that itineraries described 
by b b . . .b can be semiclassically significant in the tun- 
neling process, since smaller imaginary parts of actions 
yield larger semiclassical amplitudes. As in the Obser- 
vation 2, the same statement as the Observation 4 holds 
when the length of the sequence b b ... b is given by 2j + 1 
for a homoclinic point Wj . 

We discuss why the sequence (|62|) is bounded. 
Fig. Et c) shows the same homoclinic trajectories as in 
Fig. Etc), an( i integrable trajectories of a Hamiltonian 
H(q,p) — p 2 /2 + ke~ lq . The integrable trajectories 
are chosen so that they satisfy H(q,p) = and connect 
two infinities on the q plane, (Reg, Img) = (+00, 0) and 
(0, +00). It can be seen that the homoclinic trajecto- 
ries are along the integrable trajectories at intermediate 
time steps. In fact, as shown in Fig. I27f b). H(q,p) for 
the homoclinic trajectories almost vanishes at these time 
steps. 

The previous Fig. Elk) implies that the potential func- 
tion at the turning points, V(qj(wj)) (« V(qj-x(wj))), 
vanishes as j — > +00. Then it is expected that even for 
sufficiently large j's, H(q,p) w when the trajectories of 
Wj's are around their turning points, and that the tra- 
jectories are approximated by the integrable trajectories 
satisfying H(q,p) = 0. 

If the above expectation is correct, the imaginary parts 
of actions gained by the trajectories of ifj's described by 
a sequence b b ... b can be evaluated by the integrable 
trajectories. For our present parameter values k = 500 
and 7 = 0.005, the magnitudes of the imaginary parts 
of actions are estimated to be less than about 900 (see 
the Appendix II). This value is an over-estimated one be- 
cause homoclinic trajectories and integrable ones do not 
agree very well around the axis Reg = Img from where 
the actions along integrable trajectories are integrated. 
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FIG. 27: (a) The imaginary actions s(f~ n (wj)) plotted for 
n — 4 and j = 1 to 100. In (a),(b), and (c), the symbolic 
sequence of Wj is given by . . . OO . fob... b OO where 
b — (1, 1, v) with v = 1, 2, 3, and the length of b b ... b is 
given by 2j. The solid lines indicate the values of the imag- 
inary parts of actions in the case that the length of b b . . .b 
is one. (b) The magnitudes of the Hamiltonian H(q,p) = 

p 2 /2 + ke ,q (solid curves) along the forward trajectories of 
Wj's for j = 100. The dotted line represents the magnitude 
of the Hamiltonian at the origin (0,0). (c) The homoclinic 
trajectories already shown in Fig. 1241 c). and integrable tra- 
jectories (dotted curves) which satisfy H(q,p) = and con- 
nect two infinities of the q plane, (Reg, Img) = (+00, 0) and 
(0, +00). The dotted line sprouting from the origin represents 
Re q = Im q. The solid lines are the guides for the eye. 
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By making use of the Observations 3 and 4, we es- 
timate the imaginary parts of actions gained by the or- 
bits of homoclinic points from the symbolic sequences as- 
signed to the points. Let w be a homoclinic point whose 
symbolic sequence is given by 



. a_ 2 a-i ■ a a\ a 2 ■ 



(66) 



where a/. = (x^, yk, v^) for k G Z. We assume that vy. is 
large if a fe ^ (0,0,0). 

In the case that the symbolic sequence in ijotJj) does 
not include a consecutive part, b b . . .b, we approximate 
qk(w) for any k € Z according to the Observation 1 by 



qk(w) w (^fc7r/7) 1/2 {%k,yk)- 



(67) 



For a fc = (0,0,0), substituting {x k ,y k ) = (0,0) into (|67jl . 
we approximate qk(w) by 0. From the relation pk-i = 
<?fc - <?fc-i, we approximate p fe _i(w) by 



Pk-i(w) ~ (vr/7) 1/2 (x fc ^ fc 1/2 - x fc _ii/ fc _i 1/2 , 



1/2 



2/fe-i^fc-i 1/2 ) 



(68) 



Then due to l|55|l . ImT(pfc_i(u;)) is approximated by 

Imr(p fc _i(u))) w (71-/7) (x fc !/ fc 1/2 - a;fe-i^fc-i 1/2 ) 

x(y fc ^ 1/2 -y fc -i^_i 1/2 ). (69) 

Since the imaginary part of action gained at each time 
step is dominated by the kinetic part, as discussed at the 
Observation 3, s(w) is estimated by the sum over the 
terms in the r.h.s. of i|69|) for k > 1. 

In the case that the symbolic sequence in (|66|l includes 
a consecutive part, b b ... 6, the imaginary part of action 
gained along the itinerary described by b b ... b is negli- 
gible compared to that along the other part of the tra- 
jectory, as discussed below the Observation 4. Based on 
this fact, we approximate the imaginary part of action for 
b b ... b by a null value. This approximation allows us to 
evaluate the imaginary part of action only by the kinetic 
part also in the case of b b ... b. It is because the r.h.s. 
of jnHJ) is null when (xk-i,yk-i,Vk-i) = {xk,Vk,Vk) so 
that the imaginary part of action associated with b b . . .b 
is evaluated as a null value. As a result, whether a con- 
secutive part b b ... b is included in a symbolic sequence 
or not, the imaginary part of action is estimated only by 
the kinetic part of action. 

Finally, for the homoclinic point w whose symbolic se- 
quence is given by (|66|l . s(w) is estimated as follows: 

+00 



s(w) w y^ImT(pfc_i(w)) 



fe=l 
+00 



71 V~V 1/2 1/2 w 1/2 1/2 \ 

-Zj?kV k ' ~ Xk-iVkUKVkVk -Ifc-iffc-i)- 
' fc=l 

(70) 



Fig. [2H1 shows the imaginary parts of actions evaluated 
from actual trajectories of homoclinic points and from 
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FIG. 28: Imaginary parts of actions evaluated from (a) 
actual trajectories of homoclinic points, and (b) symbolic 
sequences assigned to the homoclinic points. In each fig- 
ure, the origin represents a null imaginary part of action 
gained by the trajectory of a fixed point associated with a 
symbolic sequence ...O O . O O O ... (O = (0, 0, 0)). The 
first to the third columns labeled by a\, a\a2, and 010203 
show respectively the imaginary parts of actions gained by 
the trajectories of the homoclinic points associated with 
...OO.OOOcnOO..., ...OO .OOO a x a-i OO..., 
and . . . O O . O O O ai 02 03 O O . . ., where a± , 02 , 03 are 
the members of {(1, 1, v) \ v = 1, 2, 3}. 



symbolic sequences assigned to the homoclinic points. 
The estimation l|70ll is based on the assumption that each 
i>k is large if ^ (0, 0, 0), however, as shown in the fig- 
ure, the estimation is still valid for small Vk- This is 
because the approximation in (|67|l are not so crude for 
small Vk , as is shown in Fig. I22f a~) . 



In each column of Fig. I28f ah the smallest s(w) is as- 
sociated with a symbolic sequence . . .0 O .0 O . . . or 
...OO.OOOb b...bO O ... with b = (1,1,1). The 
Appendix II shows such types of symbolic sequences, in- 
cluding the case that b = (—1,-1,1), attain the small- 
est s(w) in the whole candidates. In the previous Sec. 
Ill Fl we evaluated the tunneling wave functions by the 
semiclassical candidate orbits described by these types 
of sequences, typical behavior of which was illustrated in 
Fig. Etc). 
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IV. CONCLUSION AND DISCUSSION 

A. Conclusion 

We have carried out complex semiclassical analysis for 
the tunneling problem of a simple scattering map which 
creates chaotic dynamics in complex domain. Although 
classical motions in real phase space are simple, Tunnel- 
ing wave functions exhibit a complicated pattern, which 
is typically observed in chaotic systems. The wave func- 
tions were reproduced semiclassically in excellent agree- 
ment with fully quantum calculations. It enables us to 
interpret the creation of the complicated pattern appear- 
ing in the tunneling regime. Complex orbits contribut- 
ing to the semiclassical wave function are embedded in 
a hierarchical structure of an initial- value set. The hi- 
erarchical structure is a reflection of the emergence of 
homoclinic tanglement in complex phase space, i.e., the 
manifestation of complex-domain chaos. On the basis of 
symbolic dynamics constructed in the complex domain, 
phase space itineraries of tunneling orbits were related 
with the amounts of imaginary parts of actions gained by 
the orbits. Incorporation of symbolic dynamics with the 
complex semiclassical method has enabled us to discuss 
quantitatively the competition among tunneling orbits, 
and has elucidated the significant role of complex-domain 
chaos in the tunneling process of non-integrable systems. 



The chaotic dynamics in phase space is created on the 
Julia set, which includes the complex homoclinic tangle- 
ment investigated here. The trajectories on this set are 
proven to be sufficient to describe tunneling in the case 
of the complex Hcnon map [Tt| . It was numerically con- 
firmed here that this statement is correct also in our case. 
Therefore, on the basis of our present study and Ref . > 
we would like to present the notion of "chaotic tunnel- 
ing" , which first appeared in Ref. y| , as the tunneling in 
the presence of the Julia set. 

In energy domain approaches 0, 0, 0, El, to our 
knowledge, the complex-domain chaos has not been used 
explicitly in semiclassical calculations. The significant 
role of the complex domain chaos played in the time do- 
main approach should have the correspondence in the 
energy-domain ones. However, the instanton concept, 
which is intrinsic to these approaches, makes it difficult 
to see such correspondence. The reason is that even when 
one takes full account of complex classical dynamics, the 
degree of freedom of the path deformation on the com- 
plex time plane often allows one to consider complicated 
classical processes in complex domain as the composition 
of real-domain chaotic processes and instanton-like ones 
0, 0- We would like to describe the tunneling phe- 
nomena in non-integrable systems in terms of the simple 
notion, chaotic tunneling. Therefore the role of complex- 
domain chaos in the energy domain approaches is desired 
to be clarified in further studies. 



B. Chaotic Tunneling 

We further discuss the role of complex-domain chaos 
played in the semiclassical description of tunneling in 
non-integrable systems. In the present study, we adopted 
a time-domain approach of the complex semiclassical 
method. This approach is concerned with the real-time 
classical propagation, and has nothing to do with the in- 
stanton process. This means that real-domain paths are 
not connected to complex-domain paths, in other words, 
both the real domain and the other domain are invariant 
under the classical dynamics. Therefore all candidate or- 
bits to describe tunneling are always exposed to complex- 
domain chaos, not to real-domain one. In this sense, it is 
natural to consider the role of the complex-domain chaos 
in our approach. 

On our semiclassical framework, initial and final quan- 
tum states are identified with classical manifolds in com- 
plex phase space. The evolution of the manifolds is in- 
volved in the stretching and folding dynamics in the com- 
plex domain. The chain-like structure which we observed 
is nothing but the structure of the section of one back- 
ward evolved manifold, cut by the other mani- 
fold I. Our result here strongly suggests that the creation 
of the chain-like structure is only due to the emergence 
of complex-domain chaos, irrespective of the existence of 
real-domain chaos and also irrespective of the types of 
tunneling, i.e., whether energy-barrier tunneling or dy- 
namical one Wll. 



C. Discussion 

Finally, we here itemize several future problems where 
are necessary to make our theory more self-contained: 

1. We have constructed a partition of phase space in 
terms of the phase part of the gradient of a potential 
function. A similar approach can be found in the context 
of the study on a dynamical system of an exponential 
map of one complex variable |28| , where the boundaries 
of a partition correspond to the contour curves of the 
phase part of the exponential function. Genericness of 
our approach should be examined in further studies. 

2. Reproducing tunneling wave functions, we did not 
enter into details of the treatment of the Stokes phe- 
nomenon. Empirically, symbolic sequences which include 
members of the form (1, —1, v) or (—1, 1, v) with v € Z 
should be excluded from the whole candidates. In par- 
ticular, according to such empirical rule, we have ex- 
cluded from the candidates those trajectories which have 
almost null imaginary parts of actions due to the can- 
cellation between the imaginary parts gained at individ- 
ual time steps. When the condition q n +k = 1n-k-i 
is satisfied for any k > with n being fixed, where 
the asterisk denotes the complex conjugate, the imag- 
inary parts of actions integrated over the whole time 
axis become null. We observed numerically that such 
condition is satisfied by the symbolic sequences of ho- 
moclinic points such that the relation between symbols, 
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(x n +k, Vn+k, Vn+k) = (X n -k-l,-y n -k-l,^n-k-l), holds 

for any k > with n being fixed. In fact, there are an 
infinite number of symbolic sequences satisfying such re- 
lation. The criterion for whether tunnneling orbits well 
approximated by the homoclinic orbits described by such 
symbolic sequences are semiclassically contributable or 
not would be beyond our intuitive expectation based on 
the amount of imaginary parts of actions |2J| . The crite- 
rion should be given only by a rigorous treatment of the 
Stokes phenomenon. The justification of our empirical 
rule mentioned above needs the consideration of the in- 
tersection problem of the Stokes curves, and we are now 
investigating this issue [2!| . 

3. In many non-integrable open systems with the con- 
dition that V'(q) — > as \q\ — ► +oo, real-domain trajec- 
tories which diverge to infinity are indifferent, i.e., have 
null Lyapunov exponents, in contrast to the case of open 
systems with polynomial potential functions. Because of 
that, in the former systems, generaic properties of com- 
plex trajectories exploring in the vicinity of real-domain 
asymptotic region are not obvious, in spite of their semi- 
classical significant role as has been seen in our present 
study. The result of the investigation of this issue will be 
reported elsewhere. 
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APPENDIX A: ORDER OF SYMBOLIC 
SEQUENCES ACCORDING TO THE AMOUNTS 
OF IMAGINARY PARTS OF ACTIONS 

In this appendix, we give an order of symbolic se- 
quences according to the amounts of imaginary parts of 
actions estimated by the formula (|2"T|) . We consider a set 
of symbolic sequences: 

£ = {ao oi . . . a n O O . . . | n > 0, au € S' for k > 0}, 

(Al) 

where 

S' = {(1,1,1/), (-1,-1,1/) | ^eN}U{(0,0,0)}. (A2) 

For any two members of S, a and a', we introduce an 
equivalent relation '~' by 

a ~ a' 3(a) = s(a'), (A3) 

where s(cr) denotes the imaginary part of action esti- 
mated by (|27|l for the symbolic sequence a. For any two 
members of [a] and [a'], the order between them 

is defined by 

[a] < [a 1 ] & 3(a) < s(a'), (A4) 



It is easily checked that 3(a) > for any a S £ and 

the equality holds if and only if a = O O O Then 

one obtains 

[O O 0...]<[a] & a^OOO.... (A5) 

For any [a] with [a] ^ [O O O . . .], the symbolic sequence 
a can be chosen so as to take a form: 

a = ao ax - ■ .a n O O . . . (n > 0), (A6) 

which satisfies 

a n 7^ O, afc_i 7^ at for < k < n. (A7) 

For the symbolic sequence a in (|A6|) and i|A7|l . the fol- 
lowing relations hold: 

[a n O 0...}< [a n _ x a„ O O ..]< ... 
...<[a s 1 ...a„00...] = [4 (A8) 

Since, for any v G N, 

s((l,l,v) O O...) 
= «((-l,-l,t/) O O...) 
< 5((-l,-l,i/+l) O O...) 
= s((l,l,i/+l) 0...), (A9) 

either of the following relations holds: 

[(1,1,1) 0...} < [cr], (AlOa) 
[(1,1,1) 0...] = [a]. (AlOb) 

Finally, from l|A5p and (|A10|I . one obtains the following 
relations: 

[0 0...]< [(1, 1, 1) O O . . .} < [a], (All) 

for [cr] with [a]^[0 O O ...], [(1, 1, 1) O O . . .}. It is not 
difficult to check that 

[(1,1,1) 0...} 
= {b b...bO O... I 6= (1,1,1) or (-1, -1, 1)}. (A12) 



APPENDIX B: IMAGINARY PARTS OF 
ACTIONS FOR INTEGRABLE TRAJECTORIES 

In this appendix, we evaluate the imaginary parts of 
actions for integrable trajectories which satisfy T(p) + 
V{q) — and connect two infinities of the q plane, 
(Rec/, Imq) = (+oo,0) and (0, +00). These trajectories 
are included in the first quadrant of the q plane as shown 
in Fig. 1271 There are symmetric counterparts of the tra- 
jectories in the other quadrants, and the application of 
the result here to them is straightforward. 

From the relation T(p) + V(q) = 0, one obtains 

p{q) = ±iV2ke-< q2/2 . (Bl) 
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Then the action S(q,q') defined by 

SW) = f [T(p)-V(q)}^, (B2) 
can be written as 

S(q,q') = ±i2V2k I* e"< q2/2 dq. (B3) 

Jq 

We denote go = (0,0), g^ = (+00, 0), and q x ~ (x,x) 
with x > 0. The action integrated along the real axis 
from go to is evaluated immediately as 

5(g ,goo) - ±i2^kh. (B4) 

Let I be one of the integrable trajectories on the q plane 
and (x',x') be the intersection point between I and the 
axis Reg = Img. Deforming the integral path, the above 
action is represented as 

S(qo,qoo) = S(q ,q x ') +S(q x >,qoo), (B5) 

where S(qo,q x ') and S(q x i,q oc ) are integrated along the 
axis Re q — Im q and the path I respectively. 
S(qo, q X ') is represented as 

± 2^kh{[~C{w) + S(w)} + i [C(w) + S(w)]} ,(B6) 



where w — y/2j/irx' and C(w), S(w) are the Fresnel's 
functions defined by 

C{w) = / coa(TTt 2 /2)dt, (B7a) 
Jo 

S(w) = / sm{TTt 2 /2)dt. (B7b) 
From (|B4|I . (|B5|) . and (|B6|> . one obtains 



ImS(q x >, qoo ) = ±2V^7^[l-CH-5(w)](B8) 



For our parameter values k — 500 and 7 = 0.005, we 
observed that any integrable trajectory connecting infini- 
ties (+00, 0) and (0, +00) on the q plane satisfies w > 1.0. 
Since 0.6 < C(w) + S(w) < 1.4 in this range of w, we 
estimate the imaginary part of action as 

llmS^.fe)! < 450. (B9) 

Since C(w) and S(w) converge to 1/2 as 10 — > +00, 
Im S(q x >, goo) vanishes in this limit. 
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